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Abstract 

We investigate the solvent mediated interactions between nanoparticles adsorbed at a liquid- 
vapor interface in comparison to the solvent mediated interactions in the bulk liquid and vapor 
phases of a Lennard-Jones solvent. Molecular dynamics simulation data for the latter are in good 
agreement with results from integral equations in the reference functional approximation and a 
simple geometric approximation. Simulation results for the solvent mediated interactions at the 
interface differ markedly from the interactions of the particles in the corresponding bulk phases. 
We find that at short interparticle distances the interactions are considerably more repulsive than 
those in either bulk phase. At long interparticle distances we find evidence for a long-ranged 
attraction. We discuss these observations in terms of interfacial interactions, namely, the three- 
phase line tension that would operate at short distances, and capillary wave interactions for longer 
interparticle distances. 
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I. INTRODUCTION 

Nanoparticles can strongly adsorb at fluid interfaces, liquid-vapor or liquid-liquid 
For nanoparticle sizes between 1 and 10^ nm, the adsorption energy estimated from simple 
thermodynamics arguments 2| varies by several orders of magnitude, from 10 to 10^ ksT. 
.t has been shown that the adsorption energy is very sensitive to the nanoparticle geometry 
3| as well as to interfacial forces such as the line tension, which can affect signiflcantly the 
stability of the nanoparticles to remain adsorbed at the interface . Moreover, 

it has been suggested that the coupling between nanoparticle geometry and external flelds 
(electric, magnetic) may cause nanoparticle orientational transitions a, [Oj. 

The stability of nanoparticles at fluid interfaces provides a route to self assemble two- 
dimensional arrays and crystals, which have numerous applications in materials science 
[lo| . These two-dimensional structures have also attracted the interest of condensed matter 
researchers, as they provide an opportunity to investigate two-dimensional phase transitions. 



gpn nnn 



either at equilibrium j^, lul, ll2| or non-equilibrium conditions [l3|,y,ll5j. 



The interactions between nanoparticles adsorbed at fluid interfaces are expected to differ 
:'rom those in the bulk. We have discussed this issue extensively in a recent review article 
l|. Nanoparticles move on a fluctuating surface that separates bulk phases (e.g. liquid 
and vapor), which markedly differ in their permittivities and densities. These differences in 
bulk properties are expected to affect the van der Waals interactions, as well as the solvent 
mediated interactions between particles adsorbed at interfaces. Additionally, fluctuation- 
induced interactions due to the soft modes of the interface capillary waves may appear, since 
particles adsorbed at fluid interfaces act as obstacles which modify the interfacial thermal 
fluctuations. This modification of the capillary wave spectrum can lead to weak long-ranged 
attractions [l6,[l3, Explicit Ising ferromagnet calculations in two and three dimensions 
have confirmed this notion, showing the existence of attractive lon g-r anged forces between 
two points pinned at an interface between two coexisting phases 19|]. In addition to the 
interface fiuctuations, single nanoparticle studies suggest that the line tension can modify 
the wetting behavior of nanoparticles at liquid- vapor interfaces . Whether these line 
tension effects are transferred to the interaction between nanoparticles is not known. The 
experimental investigation of this question is not easy, given the difficulties in measuring line 
tension effects at the nanometer scale. Experiments of nanoparticle monolayers could help 
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to tackle these questions and could give important clues on the nature of the interactions 
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operating between nanoparticle adsorbed at fluid interfaces [2( 

In order to advance in the understanding of the nature of the interactions operating 
between nanoparticles at fluid interfaces, we have performed computer simulation and the- 
oretical investigations of the solvent mediated interactions between nanoparticles adsorbed 
at liquid-vapor interfaces. We have modeled the nanoparticles and the solvent as a simple- 
liquid mixture interacting with Lennard- Jones (LJ) type potentials. The difference be- 
tween the solvent-mediated interactions in the bulk phases (liquid and vapour) and those 
at the interface is of particular interest to us. Since solvent particles of LJ type possess a 
steep, repulsive core, the solvent-mediated interactions between nanoparticles with a like- 
wise strongly repulsive core should feature the well-known depletion interactions when their 
surface-to-surface separation is less than a solvent diameter 23|]. These depletion inter- 



actions play a very important role in determining the phase behavior of polymer-colloid 



mixtures for instance 
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261]. We expect they may be equally significant in tuning 



the interactions between nanoparticles adsorbed at fluid interfaces. Unlike the latter, the 
depletion interaction of hard-sphere models of nanoparticles and solvent, is very well under- 



stood using geometric concepts 27|] that extend the well-known Derjaguin approximation 
to a more appropriate "colloidal" limit for the depletion force ^8|. For more general solvent 
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3^. These 



models of simple liquids a number of studies exist, see e.g. Refs. 
results, however, are not completely understood theoretically. For larger surface-to-surface 
distances, i.e. away from the depletion regime, the solvent-mediated interactions between 
nanoparticles in bulk solvent exhibit oscillations that decay exponentially. These oscilla- 
tions are connected to packing effects of the solvent particles, and consequently are more 
pronounced in the bulk liquid that in in the vapor phase. 

In this work we show that the effective interactions between nanoparticles adsorbed at 
a vapor-liquid interface differ significantly from the interactions in the corresponding bulk 
phases. These interactions cannot be described in terms of a simple average of the two 
contributions (liquid and vapor). This finding applies to both the short-range depletion 
regime, and the regime of intermediate distances (surface-to-surface distances of up to 
approximately 10 solvent particle diameters). We will give a tentative interpretation of 
these differences as due to purely interfacial degrees of freedom, namely, the appearance of 
a three-phase line tension force and interfacial fluctuations. 
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The paper is structured as follows: In Sec. [TTlwe introduce the model and briefly describe 
our theoretical methods (computer simulation and reference-functional integral equations). 
In Sec. mil we present our results for the solvent -mediated interaction between the nanopar- 
ticles in the bulk phases and at the interface. The bulk results will be interpreted using the 
geometrically motivated extension of the Derjaguin approximation whereas for the interpre- 
tation of the interface results we introduce the concepts of the three-phase line tension and 
of fluctuation-induced contributions to the solvent-mediated interactions. Sec. IIVI contains 
conclusions and an outlook. 



II. MODEL AND METHODS 



A. Computer Simulations 



Molecular Dynamics simulations in the canonical ensemble were performed for a nanopar- 
ticle pair at an explicit liquid-vapor interface. We have employed a set up similar to that 
used by us in previous work 5|, |6|, |7|] . The solvent is modeled as a Lennard- Jones fluid, 



ULAe,a,ro,r] = 4e 



cr 
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cr 



UJr) 



r — Tq J \^ — ^0, 
Uh3{ess, CTs, 0, r) - ULj(e,s, a,, 0, 2.5as) (r < 2.50-, 



(1) 



(2) 



(r > 2.5a,) 

where we use the interaction strength e,, and the solvent diameter a, to deflne dimensionless 
units for the temperature T* = ksT/egs, density p* = pa^, pressure p* = pal/ess and surface 
tension 7* = 7a^/es,. The interactions between the nanoparticle and the solvent are deflned 
in terms of the shifted Lennard- Jones potential in equation ([1]) j^, 



Uns{r) 



UlA^us, CTs, CTns " CT,, r) - MLj(ens, O^s, CTns " cT,, 3.5a,) (r < a„, + 2.5a,) 



(3) 







(r > a„, + 2.5a, 



where a^, = (a„ + as)/2. We have chosen a„ = 7a, for the nanoparticle diameter. Consider- 
ing a solvent diameter of about 0.3 nm, this corresponds to a nanoparticle of ~ 2 nm, which 



is close to typical sizes of small metal passivated nanoparticles jssl . Is^ . The nanoparticle- 
solvent interaction was set to e„s = l-Se^s, giving a nanoparticle-interface contact angle 
of 95±5 degrees. Because we want to isolate the solvent contributions to the interparticle 
interactions, the nanoparticle-nanoparticle interactions Unn{r) were set to zero for r > cr^, 
i.e. beyond a hard core. All the simulations were performed for a single thermodynamic 
state, T* = 0.80 ~ 0.86Tc, where Tc is the critical temperature of the spherically truncated 



and shifted Lennard- Jones model considered here [35|]. We obtained the following coexis- 
tence densities and surface tension for this thermodynamic state: p;* = 0.732, p* = 0.019, 
liquid-vapour surface tension 7* = 0.39, and pressure p* = 0.014. 

The nanoparticles were placed on the surface of a liquid slab consisting of 3x10'^ solvent 
atoms and dimensions {x,y,z} = {41, 41, 20.5}crs. The interface plane was normal to the 
z axis. The nanoparticles did not move during the simulations. A typical computation of 
the force for a specific distance involved long runs, about 5 x 10^ time steps, with a reduced 
time step of 6t* = 0.005, which corresponds in real units to times of the order of tens of 
nanoseconds. 

The depletion force experienced by the nanoparticles at surface-to-surface separation d, 
is computed through, 

f{d) = 1 (r,, ■ (f,, - f,,)) (4) 

where Tab = (fa — i"b)/|ra — r^l is the unit vector along the axis joining the centre of mass 
of the two particles, a and b, and fj^ is the force between the nanoparticle i and the solvent 
molecules. The brackets represent an ensemble average. Similar computations were per- 
formed to obtain the potential of mean force in the bulk phases, vapor and liquid. In these 
cases the density of the solvent was adjusted to ensure that the density of the fluid far from 
the nanoparticles corresponded to the coexistence densities reported above. 

In addition to the nanoparticle pair studies, we performed simulations of nanoparticle 
arrays at the same thermodynamic conditions. We considered a low nanoparticle surface 
concentration corresponding to an area per nanoparticle of A/lNnC^) = 240, where Nn and 
A are the number of nanoparticles and the interface area respectively. Here, we replaced the 
nanoparticle hard core potential by a short range repulsive contribution. 
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Unn = 4e„, (^^\ (5) 
\r - tqJ 

where = a„ — (Js- The simulations involved typically 20 nanoparticles and 0.5 — 1 x 
10^ solvent atoms. We used long trajectories, up to 10 ns, to compute the nanoparticle- 
nanoparticle pair correlation function, g{r)^ which was later inverted to extract the total 
potential of mean force, Wt, 

WT{d) = W{d) + UUd) = -kBT\ng{d) (6) 

The solvent contribution, W{d), can be extracted by subtracting the direct nanoparticle- 
nanoparticle interaction (equation ([5])) from the total potential of mean force. In these 
simulations the nanoparticles where free to fluctuate at the interface. 

We also calculated the solvent-mediated force between the nanoparticles immersed in 
the bulk vapor and liquid state, respectively, as well as for the two additional temperatures 
T* = 1.0 {p* = 0.80) and T* = 1.25 {p* = 1.71) on the isochore p* = 0.732 (i.e. away 
from coexistence). It is quite instructive to monitor the behaviour of the solvent-mediated 
force in the depletion region upon approaching coexistence (see below). Since the solvent- 
mediated force in the depletion region will be interpreted below in the "colloidal" limit 
(large nanoparticle radius), we need the nanoparticle-solvent surface tension 7„<j, which can 
be obtained from the solvation free energy Fnp of the nanoparticle by 

471 

4rf7„, = F„p-— (7) 

This definition of the surface tension is used within the integral equation approach. In order 
to compare the accuracy of the integral equation approach with simulation data, a different 
quantity %s is evaluated which is easier to simulate, 

87rRAR%s = AFnp - AttR^ARp , (8) 

with AR chosen to be small. In the limit AR 0, the relation between both surface 
tensions is given by %s = 7ns + {,R/'^)d'^nsl dR. Note that the values for and 7„s in 
Table I have been obtained using R = cr„s for the defining surface. The difference between 
'-fns and 7„s is small. The differences between 7^5 and jns should become even smaller as 
the nanoparticle radius increases. This statement can be made quantitative by considering 
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a simple cubic form in R for the free energy of insertion, 

Fnp = 47r (^p^ + -iooR^ + kR + R^ . (9) 

Such a form corresponds to the well-known scaled particle picture, with 700 representing the 
surface tension of a particle in the limit of infinite radius, i.e., a wall, and with k, k being 
correction terms. In the recently developed idea of morphological thermodynamics [sg!] , the 
coefficients k and R correspond to the free energy coefficients of the integrated mean and 
Gaussian curvature of the nanoparticle, and no more coefficients besides p, 700, /t, R are 
needed to describe the nanoparticle solvation free energy. The derivative of equation ([9]) is 
given by, 

. r.2 27 o 



thus 7„s = 7oo + i^/{'2R), which is formally Tolman's equation, showing that the term 
{R/2)d'~fns/dR above decreases as —k,/{2R) — R/R^. By dimensional arguments, one would 
expect that \Kas/ess\ ^ 1 and \Ra1/ess\ ^ 1, thus for larger nanoparticle radii, i?, we expect 
that 7„s and 7^^ quickly converge. We note that in the limit of very small particles, with 
diameters a few times the solvent diameter, the main factor influencing the numerical values 
for the surface tension is the choice of the defining surface (see ref. |6| for an illustration of 
this effect). Although there is some arbitrariness in this choice, it is physically reasonable to 
locate the surface close to the repulsive core of the nanoparticle-solvent interaction f/„s, i.e. 
approximately at the surface of an "exclusion" sphere around the nanoparticle. Note that 
for our choice R = ans we have Uns{i^ = -R) = 0, and for smaller values of r the potential 
quickly rises. 



B. Integral Equations 

The thermodynamic properties of the solvent and the solvent-mediated interaction be- 
tween the nanoparticles are obtained through the pair correlation functions gij{r) = hij{r) + l 
in the mixture of solvent particles at density pi = p with the nanoparticles at infinite dilution, 
P2 — > 0. Here, species indices i,j are 's' for the solvent and 'n' for the nanoparticles. The 
first set of relations determining the correlation functions is given by the Ornstein-Zernike 



7 



equations. At infinite dilution, they take the following form: 



hss{r) - Css{r) = phss * Cssir) , 
hsn{r) - Csn{r) = phsn * Css{r) , 



(11) 

(12) 
(13) 



with h * c(r) = J dr'h{r)c{\r' — r|) denoting the convolution product. The second set, given 
by the general closure relations 



In gij{r) + (3Uij{r) = hij{r) - c^j{r) - bij{r) 



(14) 



requires knowledge of the bridge 
Functional Approximation (RFA) 



unctions bij {(3 = l/^ksT)). We employ the Reference 



371, l38| where accurate density functionals J-'^'^^[p{r)] of a 



hard-sphere reference system are employed to determine the bridge functions: 



Sp{r) 



p{T^)=Pj9ij{'r) 



/3^^'-^[p(r)] = /5(^-f[p(r)]-^(2),.ef[^(^)])^ 
/5^(2)'-^[p(r)]) = /3F---^(pi,p2)-$^ / drcf)''-^V;POAp.(r) 

Wf c/rdr'c;f-^(r,r';p.)Ap,(r)Ap,(r'). 



(15) 

(16) 
(17) 



Here, pj-'^'^^'^'^^ is the reference system density functional Taylor-expanded to second or- 
der in the density deviations Api(r) = Pi(r) — pi around the bulk densities pi and p2- 
Consequently, the first Taylor coefficient, the direct correlation function of first order, is a 
constant and given by the excess chemical potential of the reference system at bulk den- 
sity pf. —cf^^'^'^\r; Pi) = (3 pf^'^^^^pi) . The second Taylor coefficient is the direct correlation 
function of second order and depends only on the difference of the two position arguments: 



.(2),ref, 



(k ~ I"' I) Pi)- As can be seen from equation (fT6|l . the bridge functions 



are generated from all Taylor coefficients higher than second order in an expansion of the 
reference free energy functional around the bulk density. The hard-sphere diameter ai of 
;he reference solvent can be determined via a minimization criterion for the bulk free energy 
381], for the reference diameter a2 for the nanoparticle we choose o"2 = crn- 

The RFA usually gives very good results for the equation of state of simple liquids, see 



Table I and refs. 



371, 



381], except for the immediate vicinity of critical points. Furthermore, 



the accuracy of fundamental measure hard-sphere functionals for highly asymmetric mix- 
tures allows also a reliable determination of nanoparticle-solvent correlations hns{r). The 
effective nanoparticle-nanoparticle correlations hnn{f) are related to the solvent mediated 
potential W by l3W{r) = — ln(/;,„„(r) + 1) for r > cr„. These correlation functions can be 



determined to good accuracy outside the depletion region (r > cr„, + ag) [39|, |40|, |4l|, inside 



the depletion region the "colloidal" limit expressions of refs. 
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28( 1 can be employed (see 



below). The nanoparticle-solvent surface tension has been evaluated according to equation 
([7]) by noting that the solvation free energy Fnp of the nanoparticle is the excess chemical 
potential in the dilute limit which within RFA is given by [s^] 

= /3/x:^'^^C(p) -pjdr gUr) &„.(r) + (3 ^^[p(r)] |^(r)=p ,„.(r) ' (18) 
/^Pr'^'^^lp) = pjdr (^hnsir) [hnsir) - c^Hr)] - c(f)(r)) . (19) 

By varying the nanoparticle radius, also the surface tension 'jns defined in equation ([8]) 
could be computed and compared to the corresponding values of the MD simulation (see 
Table I ). The agreement between simulated and integral equation results is satisfactory 
and underlines the usefulness of the RFA approach for very asymmetric mixtures. We 
verified that the cubic form for Fnp (equation (Q) constitutes an excellent fit to the RFA 
numerical results for nanoparticle radii between 0.5 and 10 as- In these fits, the pressure p 
which governs the volume contribution to Fnp turns out to be almost identical to the virial 
pressure of the solvent reported in Table I . This demonstrates the good thermodynamic 
consistency of RFA. 



III. RESULTS 



A. Mean force of nanoparticles in bulk phases 

Figure [1] shows the mean force / between the nanoparticles as a function of the surface- 
surface separation d = r — an for the bulk phases, liquid and vapor, at temperature T* = 0.80 
(coexistence). The MD results for the force in the liquid phase exhibit the characteristic 
oscillatory behavior expected in a dense fluid, with regular peaks at cr^ intervals, and the 
force is attractive in the depletion region d < as- The mean force in the vapor phase 
is less structured, as expected, featuring a weak repulsive force in the depletion region, 
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indicating a preference for the nanoparticles to be solvated in the vapor phase. The mean 
force is determined in integral equation theory via /5/ = —dln{hnn{r) + l)/dr. We find 
that the force between the nanoparticles in the vapor are very well reproduced by the RFA 
integral equation in the whole interval of distances. This is consistent with previous integral 
equation studies, which showed that the hypernetted chain (HNC) and Percus-Yevick (PY) 
integral equations very accurately reproduce the structural and thermodynamics properties 



of nanoparticles in bulk vapor 42|. On the other hand those studies revealed the limitations 
of these integral equations to predict the properties in the liquid phase. The RFA integral 
equation theory employed in this work shows a dramatic improvement with respect to those 
theories. We find that the RFA theory reproduced quantitatively the simulation results for 



d> as (consistent with previous studies |39|, |40|, |4l|), but fails to capture the magnitude of 



the attractive force in the depletion region d < as- 

The mean force between the nanoparticles in the depletion region is governed by the den- 
sity distribution of the solvent spheres in the annular wedge formed between the nanopar- 
ticles. With the radius of the nanoparticle becoming large, the annular wedge induces a 
quasi-two dimensional confinement for the solvent spheres. Bulk integral equations (as em- 
ployed here) do not capture this regime well. Explicit density functional studies for hard 



spheres and large nanoparticle radii 28(] suggest the validity of the following expression for 



the depletion force in the "colloidal" limit (large nanoparticle radius): 



X « _p,_2^_,|^_L (,.2fl-.„-d). (20) 

Here, R is an effective radius of the nanoparticle (i.e. an equipotential surface). For hard 
spheres, the obvious choice would he R = ans, the radius of the exclusion sphere around 
the nanoparticle. The above expression for the depletion force follows only from geometric 



considerations on the overlap region of the two exclusion spheres [27|, |28|]. Let V{r), A{r) and 
C{r) denote the volume, surface area and integrated mean curvature of this overlap region 
(which obviously depend on the distance r between the centers of the exclusion spheres). 



According to morphological thermodynamics [36|] which we already briefly introduced in 
the discussion of the solvation free energy of one nanoparticle (equation iQ), the depletion 
potential takes the form W{r) = —pV{r) — '~fA{r) — KC{r), where p is the pressure in the 
solvent, 7 = 'jns is the nanoparticle-solvent surface tension and n is the mean curvature 
coefficient. The mean force follows as / = —dW/dr. Note that with k = 0, equation 
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(120|) turns into the well-known Derjaguin approximation for the mean force which is widely 
employed in colloidal physics. The coefficient n can actually be interpreted as a line tension 
associated with the edge of the annular wedge between the nanoparticles 43|. Obviously, 
this line tension becomes important either when x ^ in equation fl2U]) or if both p and 7 
are numerically small. 

We have tested this form for the mean force on the isochore p* = 0.732 for the tempera- 
tures T* = 0.80 (coexistence), 1.0 and 1.25 (see Table I ). In applying equation fl20|) . we have 
used R = ans = ^cFs as one would do for hard spheres and we employed the surface tension 
7 = 7ns from the RFA at the respective state points (see Table I ). For the coefficient k, we 
used the value k* = nag/ ess = —0.03 which has been determined by a fit to the simulation 
results for the depletion force at T* = 0.80 and which we use also for the other temper- 
atures.^ In Figure [2] the MD data (symbols) are compared to the theoretical expressions 
(full lines) and the Derjaguin approximation (k = 0, dashed lines). It can be seen that 
the Derjaguin approximation is almost on top of the data points for T* = 1.25, the state 
point deep in the liquid but it fails to catch the upturn of the force for d —>■ as for the two 
lower temperatures. The line tension term indeed can account for this behaviour. It is even 
the dominant term for the state point right at coexistence where both the reduced pressure 
p* ^ 0.044 (RFA) and the reduced surface tension 7*^ = 0.031 (RFA) are numerically small, 
as anticipated before. 



B. Mean force of nanoparticles at the liquid-vapor interface 

Figure [3] represents the main result of this paper. We report the solvation force between 
nanoparticles adsorbed at the liquid-vapor interface. We recall that the solvent-nanoparticle 
interactions used in our simulations result in a nanoparticle contact angle slightly larger than 
90 degrees, showing there is not a strong preference for the nanoparticles to be immersed 
in the vapor or liquid phases. Comparison of the force between the nanoparticles at the 
interface (circles in Figure [3]) with the vapor (diamonds in Figure [3]) and liquid (squares in 

^ In principle, k can be detemined from the insertion free energy of a single nanoparticle with varying radius 
and thus should be equivalent to the K-coefficient appearing in equation ([9]). However, it is not clear how 
the nanoparticle-solvent potential must be chosen for different R such that the coefficient k is unequiv- 
ocally determined. This is yet an open question in the application of morphological thermodynamics to 
soft potentials. 
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Figure [3]) counterparts shows that there is not a trivial relationship between the forces at 
interfaces and in the bulk. This is particularly evident in the depletion region, d < as, where 
the interaction between the nanoparticles at the interface is considerably more repulsive than 
in the bulk phases. The repulsive force, /* ~ 2, amounts to about 8 pN, considering the 
standard Lennard- Jones parameters for Argon (e = 0.9962 kJ/mol, a = 3.405A). Therefore, 
the simple geometric description leading to equation (!20|) for the force is not valid for the 
interface case, since it would lead to a weighted superposition of the liquid and vapour 
forces and thus to an almost zero force in the depletion region. The main conclusion from 
the data represented in Figure [3] is that the interface adds a contribution to the force that 
is not present in bulk. One possible origin for this force could be again a line tension. For 
nanoparticle separations corresponding to the depletion region, d < as, there exists a three 
phase line surrounding the nanoparticle pair. The three phase line will tend to contract 
for positive line tensions and expand for negative ones, adding an attractive or repulsive 
contribution respectively, to the total force between the nanoparticles. A repulsive force 
like the one we observe in our simulations would indicate a negative line tension. We have 
made an attempt to estimate the order of magnitude of the line tension needed to generate 
a repulsive force of the order of the one observed in the computer simulations. To this end 
we consider a phenomenological model, whereby the line tension contribution to the free 
energy is given by, = tL, where L is the length of the three phase line surrounding the 
nanoparticle pair (cf. Figure HI). Hence, within this model the force due to the line tension, 
assuming a nanoparticle contact angle of 90 degrees is given by: 



(x = 2R — an — d) . 

This contribution is very similar to the line tension contribution in equation (|20l) for the 
depletion force in the bulk. R is the radius defining the three phase line (see Figure H]). A 
natural choice for this radius would he R = a^s, as before, but one could equally well choose 
a somewhat larger value (e.g., where the nanoparticle-solvent potential Uns is minimal). 
Because at this molecular scale it is difficult to split up the force contributions due to part 
of the nanoparticle being in the vapor phase and other part in the liquid phase, we can 
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only attempt to estimate the sign and order of magnitude of the hne tension. Considering 
equation fl2T]) . the magnitude of the force when the nanoparticles are in contact {d = 0), 
and R = cr„s, we estimate a hne tension of the order of r* = —0.5. The order of magnitude 
and sign of the hne tension obtained from this analysis is similar to the values reported from 
computations of nanoparticles adsorbed at Lennard- Jones liquid-vapor interfaces, where 
line tensions of the order of r* ^ —0.3 were found js, O]. 



To further test the existence of the extra repulsive contribution between the nanoparticles 
in the depletion region, we have performed simulations of nanoparticle arrays adsorbed at 
the liquid- vapor interface (see section Hi Al for details). Figure [5] shows the corresponding 
potentials of mean force (PMF), which were obtained from integration of the forces for 
the nanoparticle pair, and from inversion of the pair correlation function (equation ([H])) 
for the nanoparticle array. The agreement between both computations in the depletion 
region is excellent, showing that the PMF obtained from the nanoparticle pair represents 
an accurate approximation for the short range interactions between particles at the low 
coverages investigated in this work. Comparison of the PMF at the interface with the bulk 
phases (c.f. Figure [5]) makes even clearer the strong impact that the interface has on the 
nanoparticle interactions. 

Finally, we discuss the long range behavior of the potential of mean force between 
nanoparticles at the liquid-vapor interface (c.f. Figure [6]). For the purpose of analyzing 
the decay of the interactions in the bulk phases we have considered the integral equation 
results. We showed above that the RFA integral equation provides an excellent approxi- 
mation to the simulation results. The potential of mean force for the bulk phases, liquid 
and vapor, decays to zero for nanoparticle separations being of the order of 6cTs, whereas 
it shows a longer-ranged behavior in the case of nanoparticles adsorbed at the interface. 
(Note that in the integration of the MD force data at the interface, we have set the potential 
to zero for the longest distance investigated in this work (lOa^)). The occurence of a long- 
ranged tail in the PMF at the interface might be connected with capillary waves. It has 
been suggested that coHoids adsorbed at fluid interfaces act as obstacles which perturb the 
fluctuations of the capillary waves through boundary conditions at the three phase contact 
line |16i]. For colloids placed at the interface at center-to-center distance r, the fluctuation 
spectrum and consequently the free energy of the capillary waves will depend on r. This 
results in a distance-dependent fluctuation force, which has been described as a variant of 
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the Casimir force (originally thought of as induced by quantum fluctuations). It has been 
shown that the fluctuation force is strongly dependent on the boundary conditions applying 
to the colloids, namely, whether the colloids are fixed or fluctuate. For fixed colloids like 
the ones investigated here and a freely fluctuating three-phase contact line, the fluctuation 
force in the limit of the capillary length A — oo is given by [3] , 

pUfiuc ^ ^ In (^1 + In + const. + O (^^^ . (22) 

This equation was derived in the context of a structureless solvent, where the spectrum of 
the interface position fluctuation is entirely governed by the capillary wave hamiltonian 44|). 
Thereby molecular effects are ignored. We expect that equation (!22|) will become a good 
approximation for long interparticle distances, since in this case the molecular nature of the 
solvent becomes secondary with respect to the fluctuations of the interface. Nonetheless, it 
is difficult to quantify the distance at which molecular effects become unimportant. This 
distance can be estimated considering the results reported in a very recent work. It has been 
shown that the dynamics of nanoscopic capillary waves on Lennard- Jones liquid surfaces 
shows very good agreement with the hydrodynamic theory down to very small wavelengths, 
of about four molecular diameters {45]. This result indicates that a continuum approximation 
may provide an accurate approach at very small length scales. 

We note that the finite size of our simulation box sets a cutoff for the capillary wave 



spectrum. Hence, equation ( 122|) can not be used directly. The 
a finite capillary length in an expansion in o"„/r is given by 17 



uU leading term in U fiuc for 



/3t/..4ln((-.n|i + l)V.„^^)+0(^), (23) 

where A ^ 1.12A. In Figure [6] we show the predicted fluctuation force for a capillary length 
X = 20 as which corresponds roughly to half the lateral box size used in the simulations and 
is thus intuitively the correct cutoff for the capillary waves. (Note that also Ufiuc has been 
set to zero at d = 10 0"^ which corresponds to r = 17 as)- We find that this force has an 
order of magnitude and decay that agrees well with the simulation results, indicating that 
capillary fluctuations could add a non negligible contribution to the interactions between 
nanoparticles at long separations. The small magnitude of Ujiuc is also due to the smallness 
of the capillary length; for a realistic capillary length of about 1 mm, the fluctuation-induced 
potential increases by a factor 4 ... 5 in the distance regime shown in Figure [61 
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IV. SUMMARY AND FINAL REMARKS 



We have reported computer simulations and integral equation studies of the interactions 
between a nanoparticle pair immersed in bulk liquid or vapor and adsorbed at the liquid- 
vapor interface of a simple liquid. We have considered a situation in which the nanoparticles 
have a contact angle of the order of 90°, hence they do not show a strong preference for the 
vapor or the liquid phase. 

The potential of mean force for the nanoparticles in the bulk phases can be explained 
well by a combination of geometric arguments for the depletion region (where the surface- 
to-surface separation of the nanoparticles is less than one solvent diameter) and of reference 
functional integral equations for larger separations. We have found that the potential of 
mean force between the nanoparticles adsorbed at the interface is significantly different 
from that obtained in the liquid and vapor phases. The interactions at the interface can- 
not be explained as a simple average of the interactions in the corresponding bulk phases. 
Specifically, in the depletion region the nanoparticles repel each other more strongly than 
in any of the bulk phases. This repulsive interaction is of the order of 2-3 ksT at zero 
nanoparticle separation. We have confirmed the existence of this repulsive interaction using 
computer simulations of low density nanoparticle monolayers adsorbed at the liquid-vapor 
interface, where the potential of mean force can be alternatively obtained from inversion of 
the nanoparticle-nanoparticle pair correlation function. 

Our analysis suggests that the line tension associated to the three phase line surrounding 
the nanoparticle pair may add a significant contribution to the nanoparticle interactions in 
the depletion region. We have estimated that the line tension needed to induce a repulsive 
force of the order of the one observed in the simulations is negative and of the order of 
^ —5 X lO^^^N. This order of magnitude and sign agrees with the line tensions estimated in 
previous investigations of nanoparticles adsorbed at liquid-vapor interfaces 5|, IS] . We have 
also shown that the interactions between the nanoparticles immersed in the bulk phases 
(vapor and liquid) decay to zero within about six solvent molecular diameters, whereas 
the interaction between nanoparticles at the liquid-vapor interface features a longer range 
attraction that appears to extend beyond ten solvent molecular diameters. The magnitude 



and decay of the attractive potential is of the order predicted from a theoretical ana! 



of the force arising from the perturbation of the interfacial capillary fluctuations 16 



ysi s 



17|. 



15 



Our work suggests that interfacial (line tensions and capillary waves) degrees of freedom 
may add a non negligible contribution to the depletion forces between nanoparticles at 
interfaces. Overall, our model of nanoparticles adsorbed at realistic liquid-vapor interfaces 
indicate that the depletion force contributes with pa 1-2 ksT to the total interaction between 
nanoparticles at interfaces. 
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MD 




RFA 




p* 


T* p* 7;, 


p* 


Ins 


7ns 


0.019 


0.80 0.014 


0.013 




-0.027 


0.732 


0.80 0.014 0.11 


0.044 


0.12 


0.031 


0.732 


1.0 0.80 -0.02 


0.83 


-0.05 


-0.14 


0.732 


1.25 1.71 -0.21 


1.75 


-0.27 


-0.35 



TABLE I: The thermodynamic states of the cut-off and shifted LJ solvent considered in this work. 
The nanoparticle-solvent surface tension ^ns is defined in equation ([7]), while the modified surface 
tension 7„s is defined in equation ([5]). "MD" refers to results from molecular dynamics simulation, 
and "RFA" stands for integral equations in the reference functional approximation. 
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FIG. 1: Solvation force between two nanoparticles immersed in bulk, vapor (diamonds) and liquid 
(squares) as a function of the nanoparticle-nanoparticle separation. The dashed lines represent 
the results from the RFA integral equation discussed in the text. Full lines are the prediction of 
the Derjaguin approximation. The error bars represent the typical uncertainty associated to the 
simulation data. 
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FIG. 2: Solvation force between nanoparticles immersed in a fluid phase of density p* = 0.732. 
Symbols represent simulation results for different temperatures, full lines represent the "colloidal" 
limit (equation (j20p . with k* = —0.03) and dashed lines are the predictions of the Derjaguin 
approximation (k* = 0). 




FIG. 3: Solvation force between nanoparticles adsorbed at the liquid-vapor interface (circles). 
The line is a guide to the eye. The solvation forces of the nanoparticles in the liquid (squares: 
simulation results and dashed line: RFA theory) and vapor (diamonds: simulation results and 
dashed line: RFA theory) phases are also shown for comparison. The error bars represent the 
typical uncertainty associated to the simulation data. 
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FIG. 4: Sketch of the geometrical construction used to compute the hne tension contribution to 
the solvation force. Rp is the nanoparticle radius, R the radius of the exclusion sphere and d the 
surface-surface separation. See text for details. 




FIG. 5: Potential of mean force for the nanoparticle pair in bulk and at the interface. Dashed lines 
represent the results from the RFA integral equation. Full lines represent the potential of mean 
force of the nanoparticle pair at the interface obtained from computer simulations. The potential 
of mean force for a nanoparticle array obtained from inversion of the nanoparticle pair correlation 
function is also shown (circles). 
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FIG. 6: The long range behavior of the potential of mean force: liquid (dashed line: RFA integral 
equation), vapor (dashed-dotted line: RFA integral equation), interface (circles and thin line: 
simulation results), and theoretical prediction from the fluctuation force equation (|23p (full line). 
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